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Abstract 

We consider a hydrodynamic system that models the Smectic-A liquid crystal flow. The 
model consists of the Navier-Stokes equation for the fluid velocity coupled with a fourth-order 
equation for the layer variable tp, endowed with periodic boundary conditions. We analyze 
the long-time behavior of the solutions within the theory of infinite-dimensional dissipative 
dynamical systems. We first prove that in 2D, the problem possesses a global attractor A in 
certain phase space. Then we establish the existence of an exponential attractor M. which 
entails that the global attractor A has finite fractal dimension. Moreover, we show that 
each trajectory converges to a single equilibrium by means of a suitable Lojasiewicz-Simon 
inequality. Corresponding results in 3D are also discussed. 

Keywords: Smectic-A liquid crystal flow, Navier-Stokes equations, global attractor. expo- 
nential attractor, convergence to equilibrium. 
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1 Introduction 

Smectic liquid crystal is in a liquid crystalline phase, which possesses not only some degree of 
orientational order like the nematic liquid crystal, but also some degree of positional order (layer 
structure). The local orientation of the liquid crystal molecules is usually denoted by a director 
field d. In the nematic state, molecules tend to align themselves along a preferred direction with 
no positional order of centers of mass. In the smectic phase, molecules organize themselves into 
layers that are nearly incompressible and of near constant width [6]. The layers are characterized 
by the iso-surfaces of a scalar function (p. A key property that distinguishes the smectic- A liquid 
crystals is that, the molecules tend to align themselves along the direction perpendicular to 
the layers. The study on the continuum theory for the smectic- A phase has a long history, see 
for instance, [UGHIISJUT]. A general nonlinear continuum theory for smectic-A liquid crystals 
applicable to situations with large deformations and non-trivial flows was established by E in [TJ. 
In [7], the following hydrodynamic system was proposed 

p t + v-Vp = 0, (1.1) 
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pv t + pv ■ Vv + Vp = V ■ (cr e + a d ), 
V • v = 0, 

V? t + v • Vip = A[V ■ (CS/ip) - KA 2 ip], 



(1.2) 
(1.3) 
(1.4) 



where 

a d = /ii(d T J D(v)d)d ® d + /U 4j D(v) + ^ 5 P(v)d ® d + d ® D(v)d), 
cj e = -f d ® d + A'V(V • d) (g) d - A(V • d)VVv2. 

In the above system, p is the density of the material, v is the flow velocity and ip denotes the layer 
variable. In the Smectic-A phase, molecule orientational direction lies normal to the layer that 
d = V</?. The scalar function p represents the pressure of the fluid, a d is the viscous (dissipative) 
stress tensor and a e is the elastic stress tensor (Ericksen tensor). As usual, -D(v) indicates the 
symmetric velocity gradient, D(v) = |(Vv + V T v). Due to the incompressibility of the fluid, 
there holds V • -D(v) = |Av. p\ > 0, \i± > and /15 > are dissipative coefficients in the stress 
tensor. The constant K > arises in the free energy (cf. [7J) and A > is elastic relaxation time. 

System (|l.l|) - ([1.4p can be viewed as the analog for the Smectic-A liquid crystal of the 
Ericksen-Leslie system [6l[9l[T9] for the nematic liquid crystal flow. Equation (jl.ip represents 
the conservation of mass, equation (|1.2p is the conservation of linear momentum, (jl.3p implies 
the incompressibility of the fluid and equation (|1.4p is the angular momentum equation. £ is 
the Lagrange multiplier corresponding to the constraint associated with the incompressibility of 
the layers such that |Vy| = 1. In order to relax this constraint, an often used approach is to 
introduce the Ginzburg-Landau penalization function /(d) = (|d| 2 — l)d (0 < e < 1) with the 
associated potential function F(d) = 4^-(|d| 2 — l) 2 such that /(d) = |j (cf. |3ll22]). Replacing 
the original Lagrange multiplier term £d in a e as well as in (|1.4p by /(d), we arrive at the 
evolution system that will be considered in the present paper: 

v t + vVv-^Av + Vp= V-(a d + a e ), (1.5) 
V •« = {), (1.6) 
P4 + v Vip = A(-AA 2 ^ + V-/(d)), (1.7) 

where 

a d = /ii(d r D(v)d)d ® d + /i 5 (D(v)d ® d + d ® D(v)d), 
a e = _/(d)®d + AV(V-d)®d- A(V-d)Vd. 

The first well-posedness result of the hydrodynamic system for Smectic-A liquid crystal flow 
mentioned above was obtained in [22J. The author considered an approximate system like ([1.5p — 
p.7p but with variable density (thus one also has a mass transport equation for p like (jl.ip ) in an 
open bounded domain il C M n , n = 2, 3. The system is subject to no-slip boundary condition for 
v and time-independent Dirichlet-Neumann boundary conditions for (p. The author derived the 
energy dissipative relation of the system and proved the existence of global weak solutions in both 
2D and 3D by using a semi-Galerkin procedure. Moreover, he described the global regularity 
of weak solutions (for large enough ^4 if n = 3) and provided a preliminary analysis on the 
stability of the system. Quite recently, system (|1.5p - (|1.7p with constant density and subject to 
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no-slip boundary condition for v but time-dependent Dirichlet-Neumann boundary data for <p 
was studied in [3|. The authors proved the existence of weak solutions that are bounded up to 
infinity time for the initial-boundary problem with arbitrary initial data. The existence of time- 
periodic weak solutions is also obtained. Assuming the viscosity ^4 is sufficiently large, the author 
studied the global in time regularity of the solution and proved the existence and uniqueness of 
regular solutions for both the initial-valued problem and the time-periodic problem. 

In our present paper, we consider the problem in the n-dimensional torus (n = 2, 3) T n := 
M n /Z n , namely, system (|1.5p - (|1.7|) is subjected to periodic boundary conditions. One of the 
possible reason for this choice is as follows. Contrary to the system for nematic liquid crystal flow 
(cf. e.g., |20j). now the equation fjl -T[) for (p is of fourth order type and thus lacks of the maximum 
principle. In particular, we lose the control of ||cL||l°°. We note that, the bound of ||c1||l°° plays 
an important role in the subsequent analysis in order to prove the regularity of solutions to 
system (|1.5p - (|1.7p (cf. Lemma 13.61 and Lemma l5.ip . Higher-order estimates of solutions can be 
obtained from some higher-order differential inequalities in the sprit of |20| . However, without 
the estimate of ||d||i,°°, we are not able to control certain higher-order nonlinear terms to derive 
the required higher-order differential inequalities. It seems that this is also necessary in order 
to complete the calculations in |22] , This difficulty can be bypassed if one additionally assume 
that the viscosity ^4 is sufficiently large (cf. Lemma 15.21 see also [3]). In the periodic boundary 
case, the key observation is that we can first obtain a uniform estimate on ||d||jj2, which by 
the embedding H 2 L°° yields the bound of ||c1||l°°- The proof relies on integration by parts, 
thus if we take the boundary conditions as in we are not able to get rid of certain extra 

boundary terms. 

The main propose the present paper is to be a first step towards the mathematical study 
of the long-time behavior of global solutions to the periodic boundary problem of system (|1.5D 
(|l,7p . In the 2D case, we are interested in the study of finite dimensional global attractors. We 
recall that a global attractor is the smallest compact attracting set of the phase space which 
is fully invariant for the dynamics and attracts all the bounded subsets of the phase space for 
large times. Thus, it is certainly a major step in the understanding of the long time dynamics 
of the given evolutional system. In particular, when the global attractor is proved to have 
finite fractal or Hausdorff dimension, then, although the phase space is infinite dimensional, 
the dynamics of the system becomes finite dimensional for large times and can be described 
with a finite numbers of parameters. This is the so called finite dimensional reduction. We 
refer to |31| for a detailed description. We will prove the finite dimensionality of the global 
attractor by showing the existence of an exponential attractor, which is a semi-invariant, compact 
set attracting exponentially fast the bounded subsets of the phase space. Moreover, it has 
finite fractal dimension and contains the global attractor. We refer to [5] and to |28j for a 
detailed introduction of this concept and for discussion on its importance. This approach has 
the advantage that, contrary to the volume contraction method (see |31|). it does not need 
any differentiability property of the semigroup. As a second step, we will study the long-time 
behavior of single trajectories, i.e., the convergence to single equilibrium. This is a nontrivial 
problem because the structure of the set of equilibria can be quite complicated and, moreover, 
may form a continuum. In particular, under our current periodic boundary conditions, one may 
expect that the dimension of the set of equilibria is at least n. This is because a shift in each 
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variable should give another steady state. Moreover, we note that for our system, every constant 
vector do with unit-length (|do| = 1) serves as an absolute minimizer of the functional E in (14.71) . 
We shall apply the Lojasiewicz-Simon approach (cf. L. Simon [30] ) to prove the convergence 
and obtain estimates on the convergence rate (see p1ll2 | [T4 ^[T5l[29l [32 l l33] and the references 
therein for applications to various evolution equations). In 3-D case, some partial results can be 
obtained. Since the L°°-estimate of d is still available, we can show the local existence of strong 
solutions for arbitrary initial data by higher-order energy estimates. Assuming the viscosity is 
sufficiently large, we also obtain the global existence of strong solution. Finally, we show that the 
global weak/strong solutions will converge to single equilibrium as in the 2D case. In particular, 
we prove the well-posedness and long-time behavior of global strong solutions when the initial 
data is close to a local minimizer of the energy E using the Lojasiewicz-Simon inequality, which 
improves the results in the literature that only the case near an absolute minimizer is considered 
(cf. [22], see also [20,32] for the nematic liquid crystal flow). 

The remaining part of the paper is organized as follows. Section 2 is devoted to some 
preliminaries and the main results of the paper. In Section 3, we prove that in the 2D case, the 
semigroup generated by our model on a suitable phase space possesses the global attractor A 
and an exponential attractor M.. This allows us to infer that A has finite fractal dimension. In 
Section 4, in the 2D case, we demonstrate that each trajectory converges to a single equilibrium 
and also find a convergence rate estimate. Finally, in Section 5, we discuss the results in 3-D 
case. 

2 Preliminaries and Main Results 

We denote the Lebesgue spaces with L p (T n ) (or simply L p ), p € [l,oo], and their norms with 
|| • 1 1 lp- When p = 2, we simply denote the L 2 -norm by || • || and its inner product by (•,•). With 
H s , s € K we indicate the Sobolev spaces H s (T n ) endowed with norm || • \\h»- To simplify the 
notations, we will denote the vector spaces (L p ) n , (H s ) n , (L p ) nxn , (H s ) nxn ... by L p and H s , 
respectively, and their norms are denoted in the same way as above. For any norm space X, we 
denote its subspace by X such that X = {w € X : fj n wdx = 0}. As customary, we introduce 
the following standard functional spaces for the Navier-Stokes equation 

H := {v G L 2 (T n ), V • v = 0}, V := {v € H^TP), V • v = 0}, V' := the dual of V. 

(•, •) denotes the duality product between V and V. The shorthand notation Dij will be used 
for the entries of the matrix D. We indicate with the same symbol C different constants. Special 
dependence will be indicated if it is necessary. Analogously, T> : K + — > M + denotes a generic 
monotone function. Throughout the paper, the Einstein summation convention will be used. 

We introduce the notions of weak/strong solutions to problem fjl . 5|) — fjl . T[> : 

Definition 2.1. (l)(v,tp) is a weak solution to problem (fT3| - (fTTT] m[0,T) xT" (T € (0,+oo)), 
i/v € L°°(0, T; H) D L 2 (0, T; V), (p <G L°°(0,T; H 2 ) fl £ 2 (0, T; H A ) and verifying 

(d t v,w) + ((v V)v,w) + y(Vv, Vw) = (a d + a e ,Vw), Vw € V, 
tp t + \r-V(p = \(-KA 2 ip + V • /(V</?)), a.e. in [0, T] x T n 
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v(0) = v , y>(0) = V9 , in T n . 

(2) A weak solution (v, (p) to problem (|1.5p ~ (|1.7p is a strong solution, if for T > 0, v € 
L°°(0,T;V) n L 2 (0,T;H 2 ), tp E L°°(0,T;H 4 ) n L 2 (0,T; H 6 ), and the system (fL5l) - (fTTl) is 
satisfied point-wisely in [0, T) x T n . 

The calculation in |22j (with different boundary conditions but the proof is the same) implies 
that system (|1.5p - (|1.7|) has a dissipative nature, in particular, the following basic energy law 
holds 

Proposition 2.1. Lei (v, (p) be a smooth solution to the system (|1.5p - (|1.7p . Define the total 
energy 

S(t) = hv(t)f + ^\\Ap(t)\\ 2 +[ F(d)(t)dx. (2.1) 
T/ien following identity holds: 

j t S(t) = -J ( W (d T J D(v)d) 2 + ^|Vv| 2 + 2 y u 5 |D(v)d| 2 ) ( i 2; 

-A||-KAV + V-/(d)|| 2 . (2.2) 

We can prove the existence of weak solutions to (|1.5p ~ (|1.7p by applying the semi-Galerkin 
approximation scheme as in |22] (cf. also [3,20,21]). The proof is similar to [3l[22] and we omit 
the details here. 

Theorem 2.1. [Existence of weak solution] Suppose n = 2,3. For any (vo, </?()) £ H x H 2 , 
system (|1.5p ~ (|1.7p admits at least one weak solution. 

A weak/strong uniqueness result was obtained in [3] for system (|1.5p - ([1.7p with different 
boundary conditions (see [22j for a statement for the system with variable density). A similar 
argument yields the same conclusion for our case: 

Theorem 2.2. [Weak/strong uniqueness] If(yi,(pi) and (v'2, ( P2) are respectively a weak and a 
strong solution of (|1.5p ~ (|1.7p in [0, T], then (vi, (fx) = (v2, ^2) almost everywhere in [0, T] x T n . 

Here are the main results of the paper: 

Theorem 2.3. Suppose n = 2. 

(1) Any weak solution to system (|1.5p ~ (|1.7p becomes strong for strictly positive times such 
that for any t > 0, 

rt+l 

Il(v,^(i)|lvx^ + J \\&v(s)\\ 2 + y(s)\\ 2 H6 ds < V(\\(v ,<p )\\ HxH 2,t), (2.3) 

T> is a positive function depending on ||vo||, H^oHi? 2 ? t an d coefficients of the system. In partic- 
ular, lim i _ i .o+ F>(t) = +00. 

(2) For any (vo,y>o) G V x H 4 , system (|1.5p - (|1.7p admits a unique strong solution. 

Theorem 2.4. Suppose n = 2. Denote the phase space % x ff 2 ; where rl = {v £ H : J T2 vcfcc = 
h} and ff 2 = {99 G if 2 : J^ 2 <pdx = c}, with h being any given constant vector in M. 2 and c is an 
arbitrary constant. 
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(1) System (|1.5|l - (|1.7p processes a global attractor A with finite fractal dimension inrlx ii 2 . 
Moreover, A is bounded inVx H and it is generated by all the complete trajectories. 

(2) System (|1.5p - (|1.7p possesses an exponential attractor M. in/H x if 2 , which is bounded in 
V x if 4 . 

Theorem 2.5. Suppose n = 2. For any vo E if, tpo E if 2 , f/ie global weak solutions to problem 
(|l,5p - (|1.7p /ias £/ie following property: 

lim (||v(t)|| H i + Mt)-¥' Q o||H4) = 0, (2.4) 

t— >-+oo 

where ip^ E if 4 is a solution to the following periodic elliptic problem: 

— KA 2 (/9oo + V • f(V(poo) = 0, x E T 2 , u>ii/i / Poodx = / podx, (2.5) 
Moreover, there exists a positive constant C depending on vo, c/?o, ty^oo A, //' s such that 

\\v{t)\\m + \\<p{t) - VooWh* <C(l + t)~™, Vt>l. (2.6) 

E (0, i) is usually called Lojasiewicz exponent and it depends on cp^. 

For any vo E V , (po E if 4 , i/ie global strong solution to problem (|1.5p ~ (|1.7p /ias i/ie same 
property (|2.4p anc? (|2.6p ZioWs /or t > 0. 



Theorem 2.6. Suppose n = 3. 

fij For any (vo,v?o) E F x if 4 , problem (|1.5p - (|1.7p admits a unique local strong solution. 

(2) For any (vo,(/?o) € V x if 4 , i/ /i4 > /^(vo,^) ^s sufficiently large (cf. (|5.1ip ). problem 
(ll.5p ~ p.7p admits a unique global strong solution. 

(3) Let (y,(p) be the weak solution to problem (jl.5p — ()1.7p on [0, +oo). T/ien i/iere is some 
T* > smc/i ffcai v E L°°(T*, oo; V) D L 2 oc (T*, °°; H2 ); V - -'° C-*, °°; # 4 ) n L 2 loc (T* , oc; if 6 ). 

^ Lei tp* E ii 2 6e a local/absolute minimizer of E(ip) (cf. (|4,7p ). For any vq E V, ipo E ii 4 
satisfying HvoHh 1 5: lj ||yo ~ ^Wh* — 1> i/iere are constants o~i,o~2 G (0)1] which depend on 
(p* and coefficients of the system such that if ||v-|| < 0"i and — V*||h 2 — °"2j ^en problem 
p.5p ~ p.7p admits a unique global strong solution. 

(5) If we further assume that jj 3 vodx = 0, then the global weak/strong solution to (jl.5[) — 
p.7p enjoys the same long-time behavior as in Theorem \2. 51 with t > 1 in (|2,6p feeing replaced 
by t >T* for the weak solution. 

Remark 2.1. Due to the periodic boundary conditions, we can easily see that the mean value of 
v and (p are conserved in the evolution: 

/ v(t)dx = / vocte, / ip{t)dx = I podx, V t > 0. 

7t™ JT n ip ip 

For £/ie safce of simplicity, by replacing v (respectively <p) with vo — J T „ vodx (respectively with 
(p — J (p dx), we shall always assume that J" Tn v dx = and jj„ ip^dx = in £/ie subsequent 
proof. Since system p.5p — (JTTTj) is invariant under a shift of <p by any constant, the transformation 
on tp will not influence all our results. However, when we shift the velocity v to make it has a 
zero mean, there will be one extra lower-order term in the equations (|1.5p and p.7p respectively. 
This difference will not influence most results we obtain except the convergence of global solutions 
to equilibria fTheorem \2.5\ and point (5) in Theorem \2. 6\) . If the mean value o/v is not zero, we 
cannot apply the Poincare inequality to obtain the decay of ||v|| H i from the convergence of ||Vv||. 
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Remark 2.2. If we simply set v = 0, system (|1.5[) — tjl.7|) is reduce to the single equation (ft = 
\(—KA 2 (p + V • /(V(/j)), which has been used to model epitaxial growth of thin films with slope 
selection in 2D, where (p denotes a scaled height function of a thin film (cf. fl7\ \23\ l ). Existence 
and uniqueness of the weak solutions as well as some preliminary results on long-time behavior 
of the solutions as time goes to infinity (like sequent convergence) was obtained in /-/&]/ . 

3 Global Attractor and Exponential Attractors in 2D 

In this section we study the long time behavior of the system (|1.5p - (|1.7p in terms of global and 
exponential attractors. As suggested by Remark 12.11 we work in the phase spaces 

$ := H x H 2 , $ x := V x H 4 

with the norms ||(v,^)||| := ||v||^ + |M|^ 2 , ||(v, <^)||| >1 := ||v||^ + |M|# 4 , respectively. It is 
obvious that <3?i is compactly embedded into <&. 

Recall the definition of the global attractor (cf. |31] ) 

Definition 3.1. Suppose X is a complete metric space. Given a semigroup S(t) : X i— y X, a 
subset Ad X is the global attractor if (i) The set A is compact in X; (ii) It is strictly invariant: 
S(t)A = A, t > 0; (Hi) For every bounded set B C X and for every neighborhood O = 0{A) of 
A in M., there exists a time T = T(0) such that S{t)B C 0{A) for all t >T. 

As far as our system is concerned, we define S(t) : $ i-> $ to be the map (vo,¥>o) l— ^ 
(y(t),(f(t)). Unfortunately, Theorem 12.11 does not guarantee that S(t) is well defined on the 
phase space $, since we are not able to prove a uniqueness result for weak solutions. We will 
refer to S(t) as a solution operator, being aware of the fact that, in principle, S(t)(vo, (fo) could 
be multi-valued due to the possible non-uniqueness. In the cases in which uniqueness holds, 
with a little abuse of notation, we will still indicate with S(t) the corresponding semigroup. As 
a consequence of the possible non- uniqueness, as it will be further explained later, we will not 
directly construct the global attractor on the phase space $ but rather on the "lifted" phase 
space of ^-trajectories. 

3.1 Dissipativity 

The following lower-order uniform estimate follows from the basic energy law: 

Lemma 3.1. Suppose n = 2,3. For vo € H, (fo € H 2 , the weak solution to (|1.5[l — (|1.7[) has the 

following uniform estimates 



v(t)\\ + \\<p(t)\\ H 2<C, t>0, 



(3.1) 



where C > is a constant depending on \\vq\\, \\(po\\ H 2, K. Moreover, 




Vv(t)|| 2 + \\-KA 2 ip(t) + V • f(d{t))\\ 2 ) dt < max 



{}A} m - (3 - 2) 



Next, we prove some dissipative estimates for the weak solutions to (|1.5p - (|1.7p . 
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Lemma 3.2. Suppose n = 2. For vq € H, ipo £ H 2 , any weak solution of ([1.5|) - (|1.7|) verifies 

||(v(t)^(t))||| < P(||(v ^ )|| $ )e- at + C, (3.3) 

where the positive constants C and a are independent on the solution and depend only on the 
coefficients of the system. 

Proof. Multiplying (jl.7p with (p and integrating over T 2 , we get 

+ ia||A^|| 2 + 4 [ \Vip\Ux = - [ (y-V)(ptpdx + ^ f \Vip\ 2 dx. (3.4) 
The righthand side of (|3.4p can be estimated as follows 

4 / \^\ 2 dx < a j |v^+4iqi ; 

-[ (V.V)^<||V||||V^|| L 4|M| L 4 < -L|| V ||2 + ||| V¥3 ||2 4 |^||2 4) 

Jj2 16\ 2 

b\ > is a small constant to be determined later. For v 6 V, we infer from the Poincare 
inequality that ||v|| < Cp ||Vv||, where the constant Cp > depends only on T 2 . For ip 6 H 2 , 
we infer from the Sobolev embedding theorem, Poincare inequality and Holder inequality that 

|T 2 r^|M| < |M| L 4 < C x \\Vip\\ < Ci|T 2 |s||V V || L 4, 
where C\ is constant depending only on T 2 . As a result, 

Al -||V^4^|||4<^ 2 |T 2 |^ [\V^dx, 

e 2 Cf\T 2 \ 2 



2 

2 



<c 2 |t 2 |||vv,|| 2 4 < A^jv^ + i^ 



Hence, we deduce that 



2 + JOIIA^H 2 + ( A - |C7 2 |T 2 ||) £ \V^dx + |Mf 
< ^HVv|| 2 + A|T 2 | + ^^A^. (3.5) 



1 (I,, „, „,„, „o ( A 5 

Cp n „ ll2 A 2 i e 2 C 1 4 |T 2 ' : 
^HVv|| 2 + ? |T 2 | + ^[ 

Multiplying (|3.5p by 52 > and adding it to the basic energy law ()2.2p . we obtain 

d 
dt 



i||v|| 2 + ^||A^|| 2 + / F(d) + ^|M| 2 1 + / [M^jA^) 2 + 2/x 5 |^d| 2 ] 

2 2 J T 2 2 J T 2 



-A A> + V • /(V^)|r + ( y " ^) l|Vv|| 2 + 5 2 KA||A^|| 2 



+5 2 (A-|c7 2 |T 2 |^ | t2 |V^x + 5 2 ||^|| 2 



(A|T 2 | + fA^) 



Take <5i , 52 that satisfying 

0' 



" 2e 2 C7 2 |T 2 |^' 2 2( ^P 4e 2 C P C 2 |T 2 [^ ' 



S 



We deduce from (13.611 that 



—II II 2 

2 im 



I||v|| 2 + ^||A^|| 2 + f F(d)dx + 
+5 2 KX\\A^ + g £ |V^<fc < 5 2 (A jT 2| + (3 . 7) 



+ ^l|Vv|| 2 
4 

,2r<4|T2|2 ~ 



Define := £(t) + ^IMI 2 - It i s eas Y to see that 

C 3 (j|v|| 2 + ||Av9|| 2 + II^H 2 + \V^dx + l\> *(t) > C 2 (||v|| 2 + l^ll^), 

where C 2 , C3 are positive constants depending on T 2 , K, e, A, 114 but not on the solution. Thus, 
we can conclude that there exist two positive constants C4, C5 depending only on T, K, e, A, /X4 
such that 

^*(t) + C 4 *(t) <C 5 . 
at 



As a result, 

l|v(t)|| 2 + \W{t)\?m < tU(*) < ^e- C ^(0) + V t > 0. 

The proof is complete. □ 



3.2 Higher-order estimates 

Next, we show that the weak solutions turn out to be regular for strictly positive times. This, will 
imply the compactness of the solution operator S(t). The following lemma plays an important 
role in the subsequent proof. It is worthwhile noting that, since the coupling in equation (|1.7p is 
weak, this result is valid both for n = 2, 3. 

Lemma 3.3. Suppose n = 2,3. We have 

^||VA99|| 2 + XK ||VAVf < C7(||Vv|| 2 + || VA^|| 2 )||VAv9|| 2 + C||Vv|| 2 + C, (3.8) 

where C is a positive constant depending on ||vo||, || Vol Iff 2 an d coefficients of the system. 

Proof. We just work in the 3-D case and it is easy to verify that the same result holds in 2D. 
Multiplying (jl.7p by A 3 ^, integrating over T 3 , due to the periodic boundary condition, we have 

^llVA^f + A^IIVAVH 2 
= - f VAV • V(v • Vip)dx - A [ VAV • V [V • /(d)] dx. (3.9) 

By the uniform estimates (|3.3|) . the Agmon inequality and Gagliardo-Nirenberg inequality in 
3D, we get 

liv^iu- < c|M|| 3 |Ml| 2 <c(||va^||3 + i), 

||VV9?|| L 3 < C(||VA^||3||A^||s + ||A^||) < C(||VA^||5 + 1). 
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Now we estimate the right-hand side of (|3.9p term by term. 



/ VAV • V(v • V(p)dx 
Jt 2 



< C|| VAVll(]|Vv|| || V^Hloo + ||v|| L fl||VV^|| L a 

< C||VAV||||Vv||(||VA^||3 + l) 

< e||VAVl| 2 + ||Vv|| 2 ||VA^|| 2 + C||Vv|| 2 . 



A 



J VAV V [V ■ /(d)] da: = ~J VAV • V [(3|V<^| 2 - l)A<p] dx 



< C||VAV||||VA(^|| +C||VAV||||V^|| 2 j oo||VA^i| 
+C||VAV||||V^|| L o ||VV^|| L 3||A ¥ p|| i a 

< e ||VA 2 ^|| 2 + C||VA^|| 2 + C||VA^|| 2 ||V^||io +C||VV^|| 2 3 ||A^||i 6 ||V^||^ 

< e||VA 2 ^|| 2 + C||VA^|| 4 + C. 

Taking e = we infer from the above estiamtes that (|3.8|) holds. The proof is complete. □ 
Denote 

Q = -KA 2 (p + V • /(d). 

By the definition of Q and the Sobolev embedding theorem, we can easily derive the the following 
estimates. 



Lemma 3.4. Suppose n = 2,3. We have ||VA<p|| < C||Q||2 + C, \\A 2 ip\\ < MQ\\ + C, 



||VAV|| < £||VQ|| +C, where C is a constant depending on \\(p\\ H 2,K,e. Moreover, ||A 3 9j|| < 
■^||AQ|| + C, where C is a constant depending on \\<p\\jj3,K,e. 



Next, we prove the following higher-order estimate for (p: 

Lemma 3.5. Suppose n = 2,3. For any v$ £ H, ifQ € H 2 , the weak solution to (|1.5|) - (|1.7|) 

satisfies 

\\<p(t)\\H* < ^pZ>(ll(vo,Po)||*), Vt > 0. (3.10) 

Moreover, if we assume in addition that ifo 6 H 3 , \\ip(t)\\jj3 can be bounded by a constant 
depending on \\vq\\ and \\(Po\\h 3 uniformly in time. 

Proof. We infer from Lemma l3.1l and Lemma [3.41 that for any r > and t > 0, 

rt+r rt+r poo 

sup/ ||VA^(r)|| 2 dr <supC / \\Q(r)\\dr+Cr < C ||Q(t)|| 2 +CV < C(l+r), (3.11) 
t>0 Jt t>0 Jt Jo 

rt+r r+oo 

sup/ ||Vv(r)|| 2 dr < / ||Vv(r)|| 2 ciT < C. (3.12) 
t>0 Jt Jo 

It follows from (j3.8|) and the uniform Gronwall lemma |31| Lemma III. 1.1] that 

||VA^ + r)|| 2 <c(l + ^\ , Vt>0, (3.13) 
which together with fj3.3j) yields (|3.10D . 
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If we assume that ipo € H 3 , then by (|3.8p . (|3.1ip . (|3.12p and the standard Gronwall inequality, 
we have 

||VA^)|| 2 

< HVA^fexp (c J (||Vv(r)|| 2 + ||VA(^(x)|| 2 + l)dr 
+C j\\\Vv(s)\\ 2 + l)exp (-C^ S (||Vv(t)|| 2 + ||VA^(r)|| 2 + l)dA ds 

< ||VA^ || 2 exp (c J (||Vv(r)|| 2 + ||VAv9(r)|| 2 + l)dr\ + C j\\\Vv(s)\\ 2 + l)ds 

< C, VtG [0,1], 

where C is a constant depending on ||vo|[, ||</?o||# 3 - Taking r = 1 in (|3. 13[) and using (|3,3p . we 
obtain the uniform estimate on ||¥?(t)||^3 for all t > 0. The proof is complete. □ 

By the Sobolov embedding theorem, we easily deduce the follow result 

Corollary 3.1. Suppose n = 2, 3. For any vq £ H, (fo € H 2 , we have 

\\V(p(t)U~ < i±lx»(||(vo^o)lk), Vt > 0. (3.14) 

Moreover, if we assume in addition that € H 3 , llV^i)!^ 00 can be bounded by a constant 
depending on ||vo|| and \\(Po\\h s uniformly in time. 

Using Corollary 13.11 we are able to derive the higher-order energy inequality in 2D. 

Lemma 3.6. Suppose n = 2. Let 

A(t) = \\Vv(t)\\ 2 +a\\Q(t)\\ 2 , 

where a > is a small constant to be chosen later (cf. (|3.24p below). We have 

^A(t) + ^||Av|| 2 + ^||AQ|| 2 <C(A 2 (t) + A(t)), V t > h > 0, (3.15) 

where t\ > is arbitrary and C is a constant depending on ||vo||, \\^Pq\\h 2 an d ti- Moreover, if 
we assume that (fo £ H s , ()3.15p holds for t > with C being dependent of \\vq\\, \\(Po\\h 3 - 

Proof. Recall the computation in pp. 1475] that V-a e = -(V-/(d))d- V F (d) + K A 2 <p\7 ip - 
KV ) ■ We note that (jl.5p can be written in the following form 

vj + vVv- yAv + VP = V • d d + (KA 2 p - V • /(d))d, (3.16) 

where P = p + V ( + F(d) ) . Using (|3.16p . we have 



--v-ll Vvll 2 = - / v t - Avdx 
2 (it 11 11 



■2 



[ (v • V)v • Avdx - — ||Av|| 2 - n\ I [V • ((d T L>(v)d)d ® d)] • Avdx 

-/i 5 / [V • ((D(v)d d + d ® D(v)d))] • Avdx 
7t 2 
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-f [(KA 2 <p - V ■ /(d))d] ■ Avdx. (3.17) 
Jt 2 

Using the periodic boundary condition and integration by parts, the right-hand side of (|3.17D 
can be manipulated as follows 

-fil / V • [(d T L>(v)d)d (g) d] • Avdx = ~m / V j(d k D kp d p didj)V ' Z V mdx 
Jt 2 J ■• 

= -jUi / V i(d k D kp d p didj)V iV jVidx = -fa / V ' i{d k D kp d p didj)S/ 'iD^dx 
Jt 2 Jt 2 

= -fix / (d k d p ViD kp ) 2 dx - 2/ii / D kp V id k d p didjV iDijdx 

JT 2 JT 2 

-2m / D kp d k d p diVidjViDijdx := -fj,i / (didjV '/ Aj) 2 dx + h + I 2 - 
Jt 2 Jt 2 

-H5 V • [(-D(v)d ®d + d® £>(v)d)] • Avdx 
Jt 2 

= ~A*5 / Vj(D ik d k dj)ViViVidx - fi 5 Vi(djD ik d k )ViViVjdx 
Jt 2 Jt 2 

= — A*5 / Vi(D ik d k dj)ViVjVidx - fi 5 / Vi{djD ik d k )ViS7iVjdx 
Jt 2 Jt 2 

= -2// 5 / Vi(D ik d k dj)ViDijdx 
Jt 2 

= -2/i 5 / (WiD ik d k ) 2 dx - 2fi 5 / Vidjd k D ik ViDijdx - 2^ 5 / djV id k D ik V iDijdx 
Jt 2 Jt 2 Jt 2 

:= -2fj, 5 / (V/ D ik d k fdx + I z + 4- 

- / [(if AV - V • /(d))d] • Avdx := / 5 . 

Jt 2 

Summing up, we have 

£4l|Vv|| 2 + ^||Avf + A*i / (d l d j V l D il ) 2 dx + 2 f i 5 [ {ViD ik d k ) 2 dx 
lot i j T 2 jj2 

5 

= / (v V)v- Avdx + V/fc. (3.18) 

By Lemma [3 .51 and Corollary 13. 1[ for any ii > 0, we have obtained the uniform estimate: 

\ H 3 + ||V^(<)||l«» < M, Vt>ti>0. (3.19) 



We now apply the Gagliardo-Nirenberg inequality, Lemma 13.41 and (j3. 19[) to estimate the right- 
hand side of (|3~T8|) . 

/ (v • V)v • Avdx < ||v|| L 4||Vv|| L 4||Av|| < ^||Av|| 2 + C||Vv|| 4 , 
Jt 2 o 

Since 

llVvfllVdll^Udll^ < C||Vv|| 2 (||AVll + 1) 
12 



< G||Vv|| 2 (||Q|| + 1) < C\\Q\\ 2 + C||Vv|| 4 + C||Vv|| 2 , 

we have 

h < [ (di4 7 -V/Dy) 2 <to + C7||Vv|| 2 ||Vd||?,co||d||| l0 o 

< t± [ (didjVtD^fdx + CWQf + CWVv^ + CWVvf. 

For after integrating by parts, we have 

h = -2/ii / DkpdkdpdiVidjViDijdx 
Jt 2 

= 2/ii / V iD kp d k d p diV idjDijdx + 4/xi / D kp V id k d p diV idjDijdx 

JT 2 JT 2 

+2/ii / D kp d k d p V idiV idjDijdx + 2/xi / D kp d k d p diV 'jV 'idjDijdx 
Jt 2 Jt 2 



ha + hb + he + hd, 



where 



/2a < T 



/ (didjViDijfdx + C|| VvfllVd^oc ||d||f,oo 
Jt 2 

< ^ / (^V^) 2 ^ + C||Q|| 2 + C||Vv|| 4 + C||Vv|| 2 , 
4 Jt 2 

/ 2fe + / 2c < C'||Vv|| 2 ||Vd|| 2 j0O ||d|| 2 j0O <C||Q|| 2 + C||Vv|| 4 + C||Vv|| 2 , 
/ M < ClIdllf.ccllAdllllVvll^^CCIIAVll'+lJIIAvllllVvll 

< ^||Av|| 2 + C||Q|| 2 + C||Vv|| 4 + C||Vv|| 2 . 

o 

As a consequence, 

h < ^ I (d i d j V l D ij ) 2 dx + ^||Av|| 2 + C7||Q|| 2 + C||Vv|| 4 + C||Vv|| 2 . 
4 Jt 2 o 

Next, 

h + h < ||Av||||Vv||||d|| LO c||Vd|| LO c ^ellAvf + CUVvfllVdll^lldll 2 ^ 
< ^||Av|| 2 + C\\Q\\ 2 + C||Vv|| 4 + C||Vv|| 2 . 

o 

h = I Qd-Avdx^llAvllllQllllV^llLoo^^HAvf + CIIQII 2 . 
Jt 2 ° 

It follows from (|3.18|) and the above estimates that 

-^l|Vv|| 2 + ^||Av|| 2 + m I (didjViDitfdx + tys I (ViD ik d k ) 2 dx 
at l J T 2 J T 2 

< CIIQf + CllVvf + CUVvll 2 , t>h. (3.20) 
On the other hand, by equation (|l,7p and integration by parts, we have 

-K I QA 2 ip t dx+ [ Q(V-f(d)) t dx 
Jt 2 Jt 2 
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-XK\\AQ\\ 2 + K [ AQ ■ A(v ■ Vip)dx 
Jt 2 

+4 I VQ • [(|V^| 2 - l)V(v • V^)]dx - 4 / VQ ■ [(IV99I 2 - l)VQ]ds 

< 7T 2 ' JI 2 

+% I VQ • [(Vy> • V(v • V<p))V<p]dx - ^ / VQ • [(Vp • VQ)V(^]dx 
e Jt 2 e Jt 2 

5 

-AJ^||AQ|| 2 + ^ J fe . (3.21) 



k=l 



The terms Ji, J5 on the right hand side of (|3.2ip can be estimated as follows. 



Ji = K I AQAv • Vipdx + 2K AQV k ViV k V ^dx + K I AQv-VA^ete 

JT 2 JT 2 Jt 2 

'■= Jla + Jib + Jlo 

where by the uniform estimate (j3.19p . Lemma [3.41 and the Sobolev embedding theorem, we get 

Jia < if||V^|M|Av||||AQ|| <^||Av|| 2 + 4aK2M2 ||Ag|| 2 , 
Jib < C||AQ||||Vv|| L 4|M| wM < C||AQ||||Vv||5||Av||5 
< ^||AQ|| 2 + ^||Av|| 2 + C||Vv|| 2 , 



Ju < ClIAQIHIvll^l^ll^M <C||AQ||||Vv||(||VA^|| L 4 + l) 
A K 

< _||AQ|| 2 + C||Vv|| 2 (||Q||+l) 



A K 

< — ||AQj| 2 + C||Q|| 2 + C|lVv|| 4 + C||Vv|| 2 . 



Next, 



J2 + J4 < C||VQ||(||V^o + l)(||Vv||||V^|| L - + ||v|| L 4 y\\ W 2,,) 

< C(||AQ||l||Q||l + ||Q||)||Vv|| 
A K 

< _||AQ|| 2 + C||Q|| 2 + C||Vv|| 2 , 

A K 

J3 + J5 < CUVQfdlV^ + l)^— ||AQ|| 2 + C||Q|| 2 . 
Inserting the above estimates into (|3.2ip . we obtain that 

4 II2WH 2 +(^K- 8aR2M2 ) ||AQ|| 2 - ^||Av|| 2 < CIIQII 2 + C||Vv|| 4 + C||Vv|| 2 . (3.22) 



dt \ [14 J 4a 

Multiplying (|3.22p by a and adding it to (|3.20p . we have 



d » / \ M4i, A ,,9 ( , T . 8aK 2 M 2 \ 
-A(t) + ^||Av|| 2 + a[XK ||AQ|| 2 



dt 4 \ A*4 

< C(l + Q)(||Q|| 2 + ||Vv|| 4 + ||Vt»|| 2 ), Vi>ti. (3.23) 

Taking 

a -A- < 324 » 

we conclude from (|3.23|) that (|3. 15[) holds. The lemma is proved. □ 
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Lemma 3.7. Suppose n = 2. For any vo £ H, ipo £ H 2 , the weak solution to (|1.5j) — fjl.7j) 

satisfies 

||(v, < C(t 2 ), V t > t 2 > 0, (3.25) 

where i 2 > is arbitrary and C(i 2 ) is a positive constant depending on ||vo||, ||^o||// 2 » ^2 anc ^ 
coefficients of the system. In particular, lim t2 _5. + C(i 2 ) = +oo. J/ we assume in addition that 
vo £ y and ifQ £ i/ 4 ; ||(v, <^)(i)||$ 1 can &e bounded by a constant depending on ||(vo,yo)||$>i 
uniformly in time. 

Proof. By Corollary 13, 11 Lemma [3,11 and the definition of a (cf. (|3.24p ). we infer that for arbitrary 
h > 0, 

r+oo 

/ A(t)dt < +oo. (3.26) 
Jti 

Since (|3.15p holds for t > t±, we can apply the uniform Gronwall lemma |31| Lemma III. 1.1] to 
get the following uniform estimate: for any r > 0, 

A(t + r) < C(ti) (l + i V Vt>ti, 

where C(ii) is a positive constant depending on ||vo||, || Vol I if 2 ) Since ti and r are arbitrary 
positive constants, we can prove the uniform estimate ()3.25p for any £ 2 > 0. 

If the initial data is more regular, namely, vo £ V and tpQ £ H A , by Lemma [3.51 we can easily 
show that ||(v, </?)(t)||$ 1 can be uniformly bounded by a constant depending on ||(vo, yo)||$i- 
The proof is complete. □ 

Corollary 3.2. Suppose n = 2. For any vo £ H, y?0 G H 2 , there exists t* > depending on 
ll v o||) IIvoIIh 2 ; suc/t that for all t > i*, i/je uiea/c solution to (jl.5|) — (|1.7() satisfies 

||(v,v)(*)lk < M, Vt>f, (3.27) 
where M is independent o/vo, Vo- 

Proof. It follows from Lemma [3.21 that there exists £3 depending on ||vo||, ||vo||.ff 2 > such that for 
all t > £3, the weak solution to (|1.5p - (|1.7p satisfies 

||(v,v)(i)||$ < Mi, Vt >t 3) (3.28) 
where M\ is independent of vo, Vo- Now, Lemma l3.1l and Lemma 13.41 imply that for t > t$, 

/•t+1 />oo 

sup / (||Vv(r)|| 2 + ||VA^(r)|| 2 )dr<C / (||Q(r)|| 2 + ||Vv(r)|| 2 ) + C < C, 
t>t 3 Jt Jt 3 



with C depending on M\. Then (j3.8|) and the uniform Gronwall lemma yield that ||(/?(i + l)||^3 < 
C, for t > i 3 . As a consequence, || Vy?(t)||L°° < M 2 for all i > i 3 + 1. For i > i 3 + 1, we fix a 
in (]3.24p with a = j^^p ■ Then (|3.15p holds with C only depending on M\ , M 2 • Applying the 
uniform Gronwall inequality once more, we have A(i) < M 3 , for t > t 3 + 2, where M 3 depends 
on Mi,M 2 . Finally, taking t* = t 3 + 2, we conclude the proof. □ 
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Proof of Theorem 12. 1L First, thanks to Lemma l3.7| we see that 



A(t) < C(t), Vt > 0, 



(3.29) 



where C(t) depends on ||(vo, Vo)!!* an d C(t) /* +oo for t \ + but remains uniformly bounded 
for t /* +oo. Integrating (|3.15p over [t,t + 1], recalling Lemma 13.11 we have 

J {\\&v(s)\\ 2 + y(s)\\ 2 H6 ) ds < C{t), Vt>0. (3.30) 
As a consequence of (|3.29p and (|3.30p . we immediately have 



On the other hand, it is easy to check that ()3.29p and f|3.30j) gives an analogous L 2 (t, t+1; L 2 (T 2 )) 
estimate of V • (a d + a e ) for any t > 0. Hence, by direct comparison with (|1.5|) — (|1.7p . we can see 
that 



The 2D smoothing property is thus proved. 

Finally, similar to |20| . we know that if the initial data are regular, the existence of a weak 
solution together with high-order estimates implies a strong solution, and by Theorem I2.2| the 
strong solution is actually unique. 

3.3 The global attractor and exponential attractors 

Lemma 13.21 and Corollary 13.21 entail that there exists a compact absorbing set in <£. If we had 
uniqueness for the weak solutions, this would be sufficient to prove the existence of the global 
attractor by using the classical theory on dynamical systems (see, e.g., |31j). We can overcome 
this difficulty essentially relying on the regularization of weak solutions to strong solutions for 
strictly positive times proved in Theorem 12.31 This implies that, for strictly positive times, we 
have enough regularity to ensure uniqueness by Theorem 12.21 As a consequence, we have the 
following weaker form of uniqueness, to which we refer as unique continuation: 

Proposition 3.1. Suppose n = 2. For any two weak solutions (yi,ip±) and (^2,^2) such that 
(vi(T), (fi(T)) = (v2(T),(f2(T)) at some T > 0, then it holds (vi,y?i) = (^2,^2) f or an U 
t > T. 

A possible way to construct the global attractor is to apply the theory of ^-trajectories in- 
troduced by Malek and Necas in |24| and later developed by Malek and Prazak in |25| (For 
other possible approaches, the reader is referred to, e.g., Ball [Tj or to Remark 13.11 in this pa- 
per). Besides, we can also use the ^-trajectory method to study the existence of an exponential 
attractor. 

For the sake of convenience, we recall some highlight points of the ^-trajectory method here. 
Roughly speaking, the ^-trajectory method consists in lifting the dynamics from the physical 
phase space to a space of trajectories with an arbitrary but fixed length £ > 0. More precisely, 
for our current problem, by ^-trajectory we mean any solution to (|1.5p - (|1.7p defined on the time 



(||v(s) • Vv(s)|| 2 + ||v(s) • V^OOH^) ds < C(t), Vt > 0. 




(||d t v( S )|| 2 + \\dMs)\\ 2 H2 ) ds < C(t), Vt. > 0. 



16 



interval [0,£]. Then, we endow the space of ^-trajectories denoted by SCt with the topology of 
L 2 (0,£;$). Note that weak solutions to (fT3 ]) -([T7r ]l lie (at least) in 

C w ([0,£];$) :={(v,<p) G L°°(0,£; : ((v, <p), (u,^))*,*/ G C([0,£]), V(u,^) G $'}, 

which makes it reasonable to talk about the point values of trajectories. 

The unique continuation property implies that from an end point of an ^-trajectory there 
starts at most one solution. Combined with the existence theorem, this implies that if (u, (ft) G Xi 
and T > £, then there exists a unique (v, if) which is a solution to (|1.5|) - (|1.7p on [0, T] such that 
(u, (/>) = (v, v)|[o,£]- Then we can define the semigroup 5^{t) on 

(y(t)(u, 0))(r) := (v(i + r), <p(t + r)), r G [0, £]. (3.31) 

From now on, without loss of generality, we will fix I = 1. Corollary 13.21 implies that there exists 
R > such that 

B x = {(v,<p) G $i : ||(v j¥ >)|k < R} C $i CC $ 

is a compact, absorbing set for the solution map S(t). Theorem 12.21 entails that the solution 
operator S(t) confined on B\ is indeed a semigroup. Let 

Bi := |J 5(t) J B 1 ' I>1 (3.32) 
t e [o,3b] 

where To > is a time such that S{t)B\ C B\ for all t > Tq and the closure is taken with 
respect to the weak topology of 3>i. Then i?i is a compact, absorbing and positive invariant set 
for S(t). Define 

m\ = {(u, (p) G Xn : (u, 0)(O) G Bt}. (3.33) 

Note that £$\ is indeed closed with respect to the topology of L 2 (0, 1; <3>). Using Corollary 13.21 
and Proposition 13.11 one can verify that all the assumptions in Theorem 2.1] are satisfied 
and as a result, the dynamical system (<5^(t), Xg) possesses the global attractor A. Next, we 
introduce the following map evaluation map 

e: L 2 (0,1;$) i-> $ defined by e((u, 0)) = (u(l), 0(1)). (3.34) 

Define B = e{ffl\). We see that B C <&i, thus 5(t) : $ — > <^ is a semigroup on B and B is 
positively invariant. If we can show that the map e is Lipschitz continuous on .%\ (which is 
indeed true, see (|3.52p below), then we can project the global attractor A back to the physical 
space <£> obtaining the usual global attractor A = e(A) for the dynamic system (5"(t),B). Since 
B is actually absorbing, A is also a global attractor in the phase space <3>. 

Remark 3.1. If one is interested only in the existence of the global attractor, one can reason 
as follows, without invoking the (.-trajectory theory. First of all, combining Theorem \2.1i \2."A 
and \2.'A we have that the restriction of the solution operator, named S(t), to the bounded sets of 
<I>1 is a semigroup. Moreover, Corollary V3.2\ give the dissipativity of S(t) with respect to the <I>i 
metric. As a consequence, the standard theory of dynamical systems gives the existence of the 
global attractor A attracting the bounded sets of <I>i but with respect to the ^-topology. Finally, 
the smoothing property implies that A is indeed the attractor for the weak solutions, since it 
attracts also the ^-bounded sets. 
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Our next step is to prove the finite dimensionality (in terms of fractal dimension) of the global 
attractor A constructed above and the existence of an exponential attractor. As anticipated in the 
introduction, the finite dimensionality of the global attractor will be deduced as a consequence of 
the existence of a finer attracting set, the exponential attractor. We recall the following (cf. [8]) 

Definition 3.1. A compact subset Ai of the phase space is called an exponential attractor 
for the semigroup S(t) if the following conditions are satisfied: (El ) The set Ai is positively 
invariant, i.e., S(t)Ai C Ai for all t > 0; (E2) The fractal dimension (see, e.g., f26l\31^ ) of Ad 
in <1? is finite; (E3) The set Ai attracts exponentially fast the image of the bounded subsets of the 
phase space <3?. Namely, there exist C, (3 > such that 

dist$(S(t)B,M) < Ce~P\ V bounded set Bc$, Vt > 0. (3.35) 

Note that, by construction, the exponential attractor, when it exists, always contains the 
global attractor. Thus, property (E2) gives that the global attractor has finite fractal dimension 
too. Besides its importance in proving the finite dimensionality of the global attractor, the 
existence of an exponential attractor is of interest in itself. In fact, it may resolve some of the 
major drawbacks of the global attractor, namely its arbitrary slow attraction, which makes the 
global attractor very sensitive to perturbation and to numerical approximation, and the difficulty 
in estimating its rate of convergence. We refer the readers to the recent survey |28j for more 
details and additional references. 

To prove the existence of an exponential attractor Ai, we first use the following existence 
theorem proposed in |10| , which gives an efficient strategy to obtain the existence of an exponen- 
tial attractor for the discrete semigroup generated by the iterations of a proper map S. Then in 
a second step, we construct the desired exponential attractor for the semigroup with continuous 
time. 

Lemma 3.8. (cf. |10| ) Let and Jtf{ be two Banach spaces such that J4f{ is compactly embedded 
into Jf?. Suppose Bj is a bounded closed subset of J4f. Let us give a map § : M\ — s> Bi such that 

Ph -SfcHjfi < L W b l -billje, V 61,62 GBi, (3.36) 

where the constant L is independent of b\ and 62. Then, the discrete semigroup {§(n), n G N} 
generated on M% by the iterations of the map S possesses an exponential attractor, i.e., there 
exists a compact set Aid C Bi such that (El) Aid is positively invariant: SAid C Aid; (E2) 
The fractal dimension of Aid * n * s finite: dim.f(Aid,J^') < M < +00; (E3) Aid attracts 
exponentially the images 0/B1 under the iterations of the map S: dist,#(S(n)Bi, Aid) < Ce~ Kn . 
Moreover, the positive constants M, C and k can be expressed explicitly in terms of the squeezing 
constant L, the size of the set M% and the entropy of the compact embedding 3tf\ CC ^ . 

In order to apply Lemma [3.81 one has to properly define the map §, together with the spaces 
ffl , and Bi. A typical choice for dissipative problems like (|1.5p - (jl.7p . would be (recall (|3.32D ) 

S := S(l), J^f := $, J#{ := $1, Bi := B\. 

Unfortunately, a closer inspection to system (|1.5p - (|1.7p reveals that the above choice is not 
completely satisfactory in the sense that proving a point-wise (in time) estimate for the difference 
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of two solutions in the norm of <&i appears to be difficult due to the highly nonlinear character 
of the problem. We overcome this difficulty by using the method of ^-trajectories to construct 
proper spaces and M\ and then verify the assumptions of Lemma 13.81 

As we did for the construction of the global attractor, we still set I = 1. Let us define 
M 1 := 3B\ (recall (pUSa]) ) and 

:= L 2 (0, 1; $), JTi := L 2 (0, 1; $i) n (W 1 ' 1 ^, 1; V) x ^(0, 1; L 2 )). (3.37) 

It follows from the Aubin-Lions compactness lemma that the embedding Jtf[ C is compact. 
We will apply Lemma 13.81 to the map S = ^(1) (see (|3.3ip ) acting on the set Bi. To this end, 
we only need to check the smoothing property (|3.36p . All the results in this subsection holds 
only in the two dimensional case. Moreover, we do not need any particular restriction on the 
values of the structural constants in the equations. Nevertheless, it would be quite interesting 
and important to find an explicit (and possibly sharp) dependence of the fractal dimension of 
the attractor with respect to the coefficients in the equations. 

Lemma 3.9. Suppose n = 2. Let (vi,^i) and (v 2 ,y> 2 ) be two solutions to problem (|1.5j) — (jl.7j) 

with initial conditions in B\. Denote v := vi — V2 and <p := <pi — <p 2 - Then, the following 
estimate holds 

ll^v|| £ i( 1)2 ;V') + \\dt<p\\L 2 (l,2;L2) + II (v,^)||l2(1,2;$i) < C|| (v, <p) |[£a(0,l;*) ■ ( 3 - 38 ) 
Proof. We know from Lemma [3.71 that for i = 1, 2 

||(vi,^)(*)lk Vt>0. (3.39) 

Then we test the equation for v with v and the equation for (p with A 2 ip, respectively. We obtain 

i|(||v|| 2 + ||A^|| 2 ) + ^||Av|| 2 + AAlAW 

= - [ (v-V)v 2 -vdx- [ (a(-ai)-Vvdx- [ {a\ - a e 2 ) -Vv dx 
JT 2 Jt 2 Jt 2 

5 

- f (vV^i + v 2 V(^)A 2 (^dx + A / V- (/(di) -f(d 2 ))A 2 £dx :=V K { . (3.40) 
Jt 2 Jt 2 ~{ 

Using estimate ()3.39p . it is not difficult to see that 



Ai 


< 


||v||^||Vv 2 || < e||Vv|| 2 + C(e)||v|| 2 , 


(3.41) 


A 4 


< 


v V(/3i l« A 2 ^ + v 2 L 4 \\V(p | L 4 A 2 <^ 






< 


e||A 2 ^|| 2 + C(6)(||v|| 2 + II^H^), 


(3.42) 


A 5 


< 


||V-(/(d 1 )-/(d 2 ))||||AV|| 






< 


C[{\\tpi\\loo +1)||A^|| + ||A^ 2 || x6 ||V^i + V^ 2 || L 6||V^|| L 6]||A 2 <^|| 






< 


e ||A 2 ^|| 2 + C(6)||^||^. 


(3.43) 



To estimate A 2 , we need to control \\af — a^W- 

af-ai = f M 1 [(djD(v 1 )d 1 )d 1 ®d 1 -(djD(v 2 )d 2 )d 2 ®d 2 ] 

19 



+/i 5 (£>(vi)di ® di - £>(v 2 )d 2 ® d 2 ) + /i 5 (di ® D(vi)di - d 2 ® D(v 2 )d 2 ) 
:= Ji + J 2 + J 3 - 

We give a detailed L 2 -estimate only for the terms J\ in the above decomposition, since for the 
other two (lower-order) terms J 2 , J3, the argument is essentially the same and actually simpler. 
We have 

Ji := (d7D(vi)d 1 )di®d 1 -(d 2 r J D(v 2 )d 2 )d 2 ®d 2 = d7 J D(v)d 1 d 1 ®di 
+ (dj - d 2 r )D(v 2 )didi di + djD(v 2 ) (didi <g> di - d 2 d 2 <g> d 2 ^ 
= Jla + Jib + Jlc- (3-44) 

The term J ia multiplied with Vv produces a nonnegative (hence negligible in the estimate) term 
since 

Jla : Vv = J la : D(y) = (d[/J(v)d 1 d 1 ® d x ) : D(v) 

= d7(-D(v))di(d 1 ®di):D(v)) = jd^/J^d!] 2 > 0. 

Then using Sobolev embedding theorem and (|3.39p . we have 

||Ji6|| < ||d7-d 2 r || LO o||D( V2 )||||d 1 (d 1 0d 1 )|| L oo <C\\(p\\ H s <C||^||| 4 ||^||| 2 , (3.45) 
\\Ju\\ < ||d 2 |[ L oo|| J D(v 2 )||(||ddi (g)di|| L oo + ||d 2 d(g)di||L<x. + ||d 1 d 2 <8> d|| L °°) 

< Cy\\ H s <C||^||| 4 ||^||| 2 . (3.46) 

As a consequence, 

K 2 < \\Vv\W\af - a d 2 \\ < e\\Vvf + e\\<pf H < + C\\(p\\ 2 H ,. 

Concerning K3, we have 

or? -erf = K[V{V ■ di) ® di - V(V • d 2 ) ® d 2 ] - A"[(V • di)Vdi - (V ■ d 2 )d 2 ] 
-[/(di) ® di - /(d 2 ) ® d 2 ] := J 4 + J 5 + J 6 . 

By the Sobolev embedding and ()3.39p . we obtain 

J 4 = ||V(V-di)®di- V(V-d 2 )®d 2 || < ||d|| H2 ||di|| L oo + ||d 2 || H 2||d|| L - 



< C\\<p\\ m <C\\<p\\* »-" 2 



1 1 

J 5 = ||(V-di)®Vdi-(V-d 2 )®Vd 2 || < ||d|| H i ||Vdi|| L «> + ||d 2 || H i ||Vd|| L <» 

< C||^|| H2 +C||^||^<C||^||| 4 ||^||| 2 , 

h < C\\?\\ H 2, 

which imply that 

K 3 < ||Vv||||^||| 4 ||^||| 2 < e||Vv|| 2 +C\m\ 2 H 2. 
Now we test the equation for <p by <p + (pt . Similar computations give 

j t m\ 2 + \\m 2 ) + iirf + wvtn* < cm 2 + \m\U- (3.47) 
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Summing (|3.40p with (j3.47j) . choosing e sufficiently small, we obtain 

^(||v|| 2 + 2||A^|| 2 + II^H 2 ) + ||Vv|| 2 + Ha^H 2 + y\\ 2 m < l7(||v|| 2 + yf H2 ). (3.48) 
By the Gronwall Lemma, for any < y — t < 2, 

l|v(y)|| 2 + \\m\\h < Ce 2C (||v(t)|| 2 + \\<f>(t)\\h) . (3.49) 
Taking t € [0, 1] and integrating (|3.48|) over [t,2], we infer from (|3.49p that 

||(v,^)(2)|| 2 +^ 2 (||(v,^(r)|| 2 i + ||«r)|| 2 )dr 

< C(||(v,^)(t)|| 2 +^ 2 ||(v,^)(r)|| 2 dr) <L7||(v,^)(t)|| 2 . (3.50) 

Integrating (|3.50p with respect to t over [0, 1], we finally obtain 

l r i 



ll(v^)lli 2 (i, 2;$1 ) + ll^||i 2(li2;L2) < / / \\(y^){r)f^ + \\d^(r)f)drdt 



o Jt 



< C||(v,$>)|| (3.51) 

It remains to estimate ||(?tv||^i( lj2; y). We use a duality argument. First, we recall that, for 
u € 2^(1, 2; V*), 



\ u \\v-(i,2-y>) = SU P 



(u, (p)dr 



where the sup is taken over the function <j) € L°°(1,2;V) such that ||0||ioon 2 ; y) = 1 an< l the 
duality pairing is between V and V. Consequently, thanks to (|3.39p and (|3.5ip . there holds 



\\dfv{r)\\ v <dr < J ||Vv(r)||dr + J (||v(r)||||Vvi(r)|| + ||Vv(r)||||v 2 (r)||)dr 
+ ^ a (l|fff(r)-^(r)|| + ||5f(r)-flfS(r)||)dr 

< cU \\(%(p)(r)\\ldrY . 
The proof is complete. □ 



Thanks to Lemma[3j51 we have verified that the map § := satisfies all of the assumptions 

of the abstract result Lemma 13.81 Therefore, the discrete semigroup {^(n), n £ N} possesses 
an exponential attractor in the trajectory space M\ endowed with the topology of J>$? = 
L 2 (0,1;<D). 

Multiplying (I3.48P with t and using the Gronwall lemma, we easily obtain that 

||(v, $0(1)11! < C f 1 \\(yMr)\\ldr, (3.52) 
J o 

which means that the map e (cf. (|3.34p for the definition) is Lipschitz continuous on Bi. This 
yields that projecting back to the phase space $ via 

M d :=e(M d ) = M d \ t=v (3.53) 
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the resulting Aid is indeed the exponential attractor for the discrete semigroup {S(n), n £ N} 
acting on B = e(Bi) (endowed with the topology of 

We note that for all the trajectories (v, if) starting from B\, there holds 

C \\d t v(s)\\ 2 + \\dMs)\\ 2 H2 ds < a 

Jo 

This, together with (|3,49p imply that the map (t, (vo,</?o)) >— > S(t)((vo,ipo)) is Lipschitz contin- 
uous on [0, 1] x B\ with respect to the R x $ metric. Thus, the desired exponential attractor 
Ai with continuous time (and on the whole phase space since B\ is absorbing) is given by the 
standard expression (see and |25j for further information) 

M := [J S(t)M d . 

te[o,i] 

The proof of Theorem 12.41 is complete. 



4 Convergence to Equilibrium in 2D 

Theorem 12.11 indicates that the total energy £(t) (cf. (|2.1j) ) is decreasing with respect to time, 
consequently, it serves as a global Lyapunov functional for system (|1.5|) - (|1.7p . The w-limit set 
of (vo, tpo) € H x H 2 is defined as follows: 

w(vo, </?()) = {(voodoo) G H x H 2 : 3 {U}^ /* +oo, such that 

(v(ti), (p(U)) -> (voo, <Poo) inHx H 2 }. (4.1) 

It follows from Lemma 13.71 and the well-known result on dynamical systems (cf. |31l Lemma 
1. 1.1]) that 

Lemma 4.1. cj(vo, <fo) is a non-empty bounded connected subset inVx H 4 . Furthermore, (i) it 
is invariant under the nonlinear semigroup S(t). (ii) £ is constant on w(vo,yo)- (Hi) w ( v O)Vo) 
consists of steady states of system (|1.5|) - (|1.7p . 

We note that the energy inequality obtained in Lemma 13.61 not only yields uniform higher- 
order energy estimates of weak solutions (cf. (|3.25|) ). but also indicates that the asymptotic limit 
points of weak solutions to problem (|1.5p — (jl.7|) actually have a special form. 

Lemma 4.2. For vo £ FL, <~Pq € H 2 , the weak solutions to (|1.5p ~ (|1.7p have the following property 

lim (||v(t)|| H i + ||Q(t)||) = 0. (4.2) 

t— >+oo 

Proof. We recall that for any t\ > 0, (|3.26|) holds. Using Lemma 13.61 and |35| Lemma 6.2.1], we 
conclude that lmi£_ s ._|_ 00 A(i) = 0, which together with Corollary 13.11 and the Poincare inequality 
leads to our conclusion. □ 

Lemma [4. 21 implies that for any initial data vo € H, ipo E H 2 , their corresponding asymptotic 
limit points (voo,^^) satisfy the following stationary problem (using the form (|3.16p ): 

-VPoo = Voo = 0, (4.3) 
-iTAVoc + V • /(V^oo) = 0, (4.4) 

ipoodx = 0. (4-5) 
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Lemma 4.3. When n = 2,3, for any cpo G L , problem 

- KA 2 p + V • f(Vcp) = 0, I tpdx= f <p dx (4.6) 

admits at least one weak solution <p, which is in fact smooth. If e is properly large, then the weak 
solution is unique. 

Proof. It is easy to verify that the energy E(cp) = ^-\\A(p\\ 2 + Jj 2 F(\7ip)dx admits at least one 
minimizer cp hi H 2 n {ip G L , f T „ ipdx = J Jn ipodx}, which is a weak solution to problem (|4.6p . 
Moreover, for any weak solution cp to (|4.6p . we have 

A||A</,|| 2 + \f \V<p\ 4 dx = \\\V<p\\ 2 < -L /" |V0| 4 <ix + -L|TT|, 

which together with the Poincare inequality ||^>|| < Cp (||V0|| + | Jj„ <j)dx\) implies that ||</>||#2 
can be bounded by a constant depending on | Jj„ ipodx\, e, Cp and |T n |. A bootstrap argument 
yields that (p is actually smooth and for m G N, ||</>||#m can be bounded by a constant depending 
on | Jj n ipodx\, e > Cp and |T n |. Finally, let (pi and (p 2 be two solutions of ()4.6|) . using the fact 
that J T n(/(V0i) — f{S/cp 2 )) • V(0i — <p 2 )dx > 0, then we infer from the Poincare inequality that 
A'||A(0i - (j) 2 )\\ 2 < ^||V(^»i - (p 2 )\\ 2 < if ||A(0i - (p 2 )\\ 2 . As a result, if e > C P A"i then 
(f>i = <p2- (We refer to [18] for a similar problem but with different boundary conditions) □ 

4.1 Convergence to equilibrium 

Lemma [4 . 2 1 yields the convergence of velocity field v. In what follows, we study the convergence 
for ip. First, e > CpK~^, we infer from Lemma 14.31 that u;(vo,c^o) consists of a single point 
(0,^00) where 9300 is the unique solution to (|4.4p - (|4.5|) . However, if e < CpA~2, we lose the 
uniqueness of steady states. Alternatively, we shall use the Lojasiewicz-Simon approach. Denote 
A = —A with D(A) = {(/> G H 2 , Jj n (pdx = 0}. Then A is self-adjoint and positive definite. Let 
Ha be the dual space of = {cp £ H , Jj n (f)dx = 0}. Then the norm on Ha is given by 

Uf A = I rn cl>A- 1 4>dx=\\A-h\\ 2 - 

We introduce the following Lojasiewicz-Simon type inequality: 

Lemma 4.4. Suppose n = 2,3. Let ip be the critical point of energy 

E(<p) = ^\\Av\\ 2 + [ F{Vip)dx. (4.7) 

/ Jjn 

Then, there exist constants (3 > 0, 6 G (0, ^) depending on ip such that for any ip G H 3 with 
\\<p — ip\\ff2 < (3 and jj„ ipdx = L n ipdx, there holds 

|| - AAV + V • f(Vip)\\ A > \E{p) - E(iP)\ l ' d . (4.8) 

Proof. Slightly modifying the arguments in |13U29| . we can easily prove that there exist constants 
Pi > 0, G (0, ^) depending on ip such that for any cp € H 3 with \\ip — ipWu'^ < @i and 
jj n ipdx = J T „ ipdx, ()4.8p holds. Next, we slightly relax the smallness condition and show that 
(|4,8p still holds if one only requires that ip falls into a certain .^-neighborhood of ip. For any 
cp G H s satisfying J* T „ ipdx = J T „ ipdx, using the regularity theory for elliptic problem, we can 
see that 

\\<P-iJ>\\h* <M||A 2 ^-^|U, (4.9) 
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where M is a constant independent of ip. On the other hand, if ||y — ip\\u 2 — 1 (which implies 
that Hvll/f 2 < HV'lliT 2 + 1)) then by Sobolev embedding theorem, we get 

||v • f(Vip) - v • /(v^)|U < dy - i>\\h*, 

where C\, C 2 depend on ||^||#2 an d WpWh 2 (by our assumption, the later one can be bounded by 
using only HV'llja" 2 )■ A s a consequence, there exists a (sufficiently small) f3 € (0,1] independent 
of if, such that if \\<p — ip\\}j2 < (3, then 

||v • f(y<p) - v • /(W)|U + \E(<p) - Em 1 - 9 < §jf ■ (4.10) 

Now for any 99 € if 3 satisfying Jj n ipdx = J" Tn tpdx and \\cp — < ft-, there are only two 

possibilities: (i) If \\<p — ip\\ij3 < then (|4.8|) holds, (ii) If ||</3 — V'll// 3 > /^l, noticing that ip 
satisfies (|4.5p . we deduce from f|4.9j) and (|4.10p that 

H-jcaV + v-zcv^iu = ||-^aV-^) + v-/(v^)-v-/(v^)|U 

> K\\A 2 (<p - ip)\\ A - ||V • f(Vip) - V • /(VV)|U 

> ^Ib - ^|Ih3 - ||v • /(v^) - v • /(W)IU 

> ^>l^)-w)[ 1 - fl . 

The proof is complete. □ 

For any initial data (vo, ^o) G H x H 2 , it follows from Lemma 13.71 that ||</?||#4 is uniformly 
bounded for t > t 2 > 0. Therefore, there is an increasing unbounded sequence {tj}j<=N an d a 
function 9300 € H A satisfying (|4,4p - (|4.5p such that 

lim \\ip(ti) - (PooWh'-i = 0, lim 6(U) = E((poo). (4.11) 

ti— »+oo tj.— > + oo 

By ()1.7p . we have 

\\<p t \\ < ||v • V(p\\ + ||Q|| < || Vv|| 11^11^2 + ||Q||. (4.12) 

We first exclude the trivial case, i.e., that there exists a to > such that £(to) = Efaoo). In this 
case, for all t > t Q , we deduce from ([22]) that ||Vv(t)|| = \\Q(t)\\ = 0. It follows from (j4~T2~P that 
for t > to, \\<Pt\\ = 0- Namely, ip is independent of time for all t>t$. Due to (|4.1ip . we conclude 
that (p(t) = ipoo for t > tQ. In this case, there is nothing else to prove. 

Therefore, without loss of generality, for all t > 0, we suppose that £{t) > E^tpoo). For 
arbitrary t > 0, we know that <p 6 L 2 (t,t + 1;H ) n -ff 1 (t,t + 1;L 2 ) which implies that 9? £ 
C([t,t + 1], i? 2 ). Due to this continuity, by a standard contradiction argument (see |15|). we can 
prove that there is a (sufficiently large) to > such that for all t > to, \\<p(t) — Poo\\h 2 < P- 
Namely, for all t > to, <p>{t) satisfies the conditions in Lemma 14.41 Apply Lemma 14.41 (|2.2p and 
the Poincare inequality, we obtain 

~{8{t) - E{ Voo )) e = -9{£{t) - E{ Voo )) d ~ l U{t) > ^ l|Vv| ' 2 + A l|Q| ' 2 



^V-W -vroo/y -v-w ~Vroo„ || v ||2(l-0) + || g|| ^ 

> C(||Vv|| + IIQII), Vt>t , (4.13) 
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which implies that J^°°(||Vv(t)|| + \\Q(t)\\)(1t < +00, and by KW . J+ 00 ||^(r)||dr < +00. 
This easily yields the convergence of ip(t) in I? as t —> +00. Since tp is compact in H s , we infer 
from (|4.1ip that limt_ 5 . +00 \\(f(t) — ^^11^3 = 0. By the Sobolev embedding theorem, we have 

AlAVW-AVooll < ||Q(t)|| + ||V-/(V^(t))-V-/(Vv o)|| < \\Q(t)\\+C\\<p(jb)-<Poo\\H*, (4-14) 

where C depends on ||<£>(i)||ij3 and ||y>oo||.H'3. As a consequence, we can conclude from (|4.2p and 
the .ff 3 -convergence of if that 

lim \\<p(t)-<p oo \\ H A = 0. (4.15) 

t— >+oo 

4.2 Convergence rate 

It remains to prove the convergence rate. This can be done in two steps: the first consists in 
obtaining, via the Lojasiewicz-Simon inequality (cf. e.g., |14j). the convergence rate for the lower 
order terms. In the second step, we will use the energy method to obtain the convergence rate 
for the higher order terms. From Lemma [4.41 and (|4.13p . we have 

j t (£(t)-E( Voo )) + C(£(t)-E( Voo )) 2 ^ <0, Vt>t , (4.16) 
and as a consequence, 

£{t) -Efaoo) < C(l + ty^, Vt> to- 
Integrating (|4.13|) on (i, +00), where t > to, it follows from (|4.12p that 

r+co r+00 g 

/ y t (r)\\dT<C (||Vv(r)|| + \\Q(t)\\)cIt < C(l + ( 4 .17) 
Jt Jt 

which implies 

\\<p(t) -V00W < c(i + ty^, t>t . (4.18) 

It follows from the basic energy law (|2,2p and (|4.5p that 

|y(*) + y l|Vv|| 2 + A||Q|| 2 <0, (4.19) 

where 

V(t) = Jl|v(t)|| 2 + ^||A v ?(t)-Av9 00 || 2 + f [F(V<p(t))-F(V<p x ) + V-f(V<Poo)(<p(t)-<Poo)]dx. 
A direct calculation yields 

[F(y<p) - F(Vv9oo) + V • f(y<Poo)(<P ~ foo)]dx 

= y^d-doo). [(idi 2 +id 0O i 2 +d-d 0O )(d-d 0O )]dx 
if 1 

+^2 y T2 ( d ~ d °°) ' [l d oo| 2 (d - doo) + doo(d + doo) • (d - doo)] dx - _||d - d c 
Thus, we have 

f [F(X7<p) -F(V^ 00 ) + V-/(V^ 00 )(^-^ 00 )]a!x < C\\ - Vp^f, 



/ 



|2 
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which together with the Poincare inequality implies 

y(t) > ^||v|| 2 + ^\\A<p - A^oof - C\\<p - ifocf- (4.20) 
On the other hand, it follows from (|3.43p and ()4.14p that 

K\\A 2 <p- AVooll < IIQII + c\\tp- v3oo||h2 < ||Q|| + C||AV- aVoo||5||v3-v?oo||5 

< ||Q|| + —\\A 2 <p - AVooll + C\\(p - ipoo\\, (4.21) 

which yields 

y(t) < i||v|| 2 + y||A99- Aipccf + CWVtp- V^ooll 2 



< c||Vv|| 2 + c||aV - aVooIIIIp - + c||aV - AVooll^lb - Vooll 2 

< C|[Vv|| 2 + C||Q|| 2 + C\\ip-cp 00 f. (4.22) 
For t > to > 0, Lemma [3.71 implies that A(t) < C that combined with (|3. 15[) yields 

^A(t) < CA(t). (4.23) 
It follows from (gig]) , (gig]) and (|4T20|) ([4T23]) that there exist constants M 1 ,M 2 > such that 

-[y(i) + M]A(t)] + M 2 [y(t) + AfiA(t)] < C\\<p(t) - <Ax>|| 2 < C(l + t)"^, V t > t . (4.24) 
By a similar argument as in [33], we conclude from (|4.24p that 

,(t)+MxA(t) < [y(to) + M 1 A(to)]e^°-') + Ce- M2t / e M2r (l + r)"i^dr 

Jto 

29 

< C(l + *)~i=a*, yt>t . (4.25) 



y( 



Finally, from ()4.18p . (|4.2ip and (|4.25p we obtain the required estimate 

\\v(t)\\ H i + \\<p(t) - (pooWm <C(l + t)-^ss, Vt>t . (4.26) 

5 Results in 3D 

Lemma 5.1. Suppose n = 3. We have 

^A(t) + ^||Av|| 2 + ^||AQ|| 2 <C,(A 3 (t) + A(t)), V t > h > 0, (5.1) 

where t\ > is arbitrary and a > 0, C* > are constants depending on \\vq\\, ||v?o||# 2 anc ^ 
Moreover, if we assume that (po 6 -ff 3 , ()5.ip fooWs /or t > u>ii/i C* 6em<? dependent of \\~vq\\, 
\\<Po\\h3- 

Proof. Using Lemma 13.51 and Corollary 13. 1[ we modify the calculations in Lemma 13.61 by using 
the 3D version of embedding theorems. It is not difficult to see that we are still able to choose 
a > sufficiently small in A(t) such that our conclusion holds true. □ 
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The existence of local strong solution for arbitrary viscosity [14 > is a direct consequence 
of Lemma 15.11 



Theorem 5.1. Suppose n = 3. For any (vo, </>o) & V x H 4 , problem (|1.5|) - (|1.7p admits a unique 
local strong solution. 

Proof. Since ipo £ i/ 4 , we have uniform estimates for ||y||//3 and ||V</?||l°°- (|5.1|) is valid for 
i > 0. Considering the ODE problem: 

±Y(t) = C*[(Y(t)f + Y(t)], y(0)=A(0), (5.2) 

we denote by 7 = [0,T max ) the interval of existence of the maximal solution Y(t). We thus have 
lim f ^ T - Y(t) = +00. It easily follows that for any t £ I, < A(i) < ^(i). Consequently, 
A(i) exists on I. This and Theorem 12.21 imply the local existence of a unique strong solution of 
problem dUSI) dUTJ) . □ 

Proposition 5.1. Suppose n = 3, vo £ V, ipo £ H 4 . For any R £ (0, 00), whenever ||Vvo|| 2 + 
||Q(0)|| 2 < R, there is a small constant £q £ (0, 1) depending only on R and coefficients of the 
system such that either (i) Problem (|1.5p - (|1.7p has a unique global strong solution (v,ip), or (ii) 
there is a T\ £ (0, +00) such that £{T\) < £(0) — Eq. 



Proof. The proof follows from the argument in [20] for simplified nematic liquid crystal model. 
A statement was also given for the Smectic-A system with variable density in |22j without proof. 
For the convenience of the readers, we sketch the proof here. We suppose that (v, ip) is a 
weak solution with initial data (vo,<^o) such that ||Vvo|| 2 + ||Q(0)|| 2 < R. Then C* in (|5.ip is 
determined by R. Moreover, due to Corollary 13.11 we have uniform estimate on ||V^(<)||l°° f° r 
all t > which only depends on R. We fix through R the constant a in the definition of A(i). 
Consider the ODE problem (|5.2p with Y(0) = max{l, a}R > A(0). Let Y(t) denote the unique 
maximal solution defined on [0,T max ). The time T max is determined by Y(0) and C* in such a 
way that it is increasing when Y(0) is decreasing. Now we take 



t = \T max (Y(0),C*), E = ^mm{^,\} 



If (ii) is not true, we have £(t) > £(0) — Eq for all t > 0. From the basic energy law (|2.2p . we 
infer that 

rto poo 

/ A(t)dt < / A(t)dt < keq, with k = max{l, q} max{2/i7 1 , A -1 }. 
Jo 

Hence, there exists at, £ [^,<o] sucri that A(t*) < ^f 2 - < Y(0). Restarting the flow (f5T2|) from 
£*, we infer from the above argument that A(t) remains bounded at least on [0, C [0, t* +^o] 
with the same bound as that on [0, to]. As a consequence, an iteration argument shows that A(t) 
is bounded for all t > 0. The proof is complete. □ 

As an immediate consequences of the above result, we can prove (1) eventual regularization 
of weak solutions and (2) the well-posedness of strong solutions near the absolute minimizers of 
energy E (cf. g^D) (cf. HUES]). 
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Corollary 5.1. Suppose n = 3. 

(1) Let (v, if) be the weak solution to problem (|1.5|l - (|1.7|l on [0, +oo). Then there is some 
T* > such that v G L°°{T*, oo; V) n L 2 oc (T*, °°; h2 )> V e £°°(T*, oo; H 4 ) n £ 2 oc (T*, °°; # 6 )- 

fSj Lei (£>* € H 2 be an absolute minimizer of E(tp) in the sense that E(ip*) < E(cp) for all 
if G H 2 . For any vq G V, ipo G H 4 satisfying HvoIIh 1 < 1 an d WfO ~~ V 9 *!!// 4 < lj i^ere is a 
constant a which may depend on coefficients of the system and ip* such that if ||vo|| < a and 
\\ifQ — (p*\\ij2 < a, then problem (|1.5|) - (|1.7|) admits a unique global strong solution. 

Next, we improve the second part of Corollary 15.11 by proving the well-posedness of strong 
solutions close to local minimizers of the energy E: 

Theorem 5.2. Suppose n = 3. Let <p* G H 2 be a local minimizer of E(tp) in the sense that 
E(<p*) < E(ip) for all ip € H 2 satisfying ||(/? — 9?*||#2 < 5. For any vo £ V, ipo £ H 4 satisfying 
II v IIh 1 < 1 an d H^o ~~ V 9 *!!// 4 — 1; there exist constants o\,oi £ (0,1] which may depend on 
coefficients of the system and ip* such that if ||vo|| < o~\ and ||<^o — c / 9 *||h 2 — °"2> ^en problem 
(ll.5D - (|1.7p admits a unique global strong solution. 

Proof. Without loss of generality, we assume 5 < 1. By Cj, i = 1,2, ... we denote constants that 
only depend on y?* and on coefficients of the system. If HvoHh 1 — 1 an d ll^o — 9 9 *||h 4 — ^ is 
not difficult to see that ||Vvo|| 2 + ||Q(0)|| 2 < R, where R depends only on ip* . Fix this R, using 
Proposition 15. 11 we can also fix the critical constant £q determined by R. It follows from Lemma 
13.11 and Lemma 13.51 that ||v(i)|| and ||¥?(£)||#3 are uniformly bounded (by a constant depending 
on <p*). Since £ is decreasing, we can see that 

0<£(0)-£(t) = I||vo|| 2 -^||v(t)|| 2 + J E(v9o)-^(i)) < ^l|v || 2 + E{<p Q ) - E(<p(t)) 

< \ho\\ 2 + CiMt)-<p \\ H 2. (5.3) 

First we require a\ < min 1^£q , lj- Let /3 denote the constant depending only on ip* provided 
by Lemma 14.41 If we are able to prove 

\\ip(t)-ip \\ H 2 <w :=mmh,^,5\ , Vi > 0, (5.4) 



then we can infer from (|5.3p that 

£{t) > 5(0) -e , Vt>0, (5.5) 

and our conclusion immediately follows from Proposition 15.11 We prove (j5.4|) by a contradiction 
argument. Assume a 2 < ^. Let t» denote the smallest and finite time for which || <£>(£*) — 
V*!!// 2 > Without loss of generality, we can assume that £(i) > E{ip*) t G [0,f*). In 
fact, if there exists t** G (0,i*) such that £(£**) = E((p*), since y>* is the local minimizer and 
||y(i**) — < / ? *||h 2 < tu < 5, we can see that v(t**) = and </?(£**) = 92**, where ip** is also a local 
minimizer (possibly different from <p*) satisfying (|4.4p - (|4.5p . Due to the uniqueness of strong 
solution, the evolution starting from i** will be stationary and hence contradicting the definition 
of i*. Thus, let £{t) > E(<p*) for t G [0,t*). We observe that the conditions in Lemma 14.41 are 
fulfilled with ip* , on the interval [0,i*). In analogy with (|4.13p . we obtain 

- j t (£(t) - E(ip*)) e > C 2 (||Vv|| + IIQII), V t G [0,Q, (5.6) 
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where C2 depends on #,/i4, A. Using (|4.12p . we have 

f Q * \Wt{t)\\dt < C 3 (£(0)-E(^)) e <C 3 (^j e \\vo\\ 2e + C 3 \E(w)-E(^)\ e 
< C4\\v \\ 2e + C 5 \\(p - (p*\f H 2. 



As a result, 



\ip{U) - <p \\jj2 < \\if(t* )- <f \\ H 2 + \\ifo - <f*\\ H 2 

< C\\ip(t*) - ^olUfsllvi**) - <£oll 3 + llvo - (p*\\m 



2 



Taking 



< C(\\<p(t*)\\ HS + \\<p \\ H 3)z U \\Mt)\\dt) +\\<Po-<P*\\h> 

< C 6 [ \\-vq\\t + \\ip -tp*\\^ 2 ) +\\ip Q -ip*\\ H 2, (5.7) 



-1 < min [\el 1, (^) " | , a 2 < min | H , (jL\ ° \ , (5. 



we easily infer from (|5.7p that — V* ll/f 2 — 1 117 < which leads to a contradiction with the 

definition of t*. Hence, we have shown that (|5,4p holds for all t > 0. The proof is complete. □ 

Remark 5.1. TTie above proof indicates that if the (regular) initial data are properly close to 
certain local minimizer, then the global strong solution will remain in the neighborhood of this 
local minimizer for all time. The conclusion is also true for the case with absolute minimizer in 
Corollaru \5.1\ If the minimizer is isolate, then we obtain the stability of it. 

Lemma 5.2. Suppose n = 3. Denote Ai(i) = ||Vv(i)|| 2 + ||Q(t)|| 2 and Ai(f) = Ai(i) + 1. For 

any /14 > 0, we have 

lAxM+^-M^lMit)) ||Av|| 2 + (xK-M 2 ^(l + ^)A 1 (t)\ \\AQ\\ 2 

< M 3 (l + Mj 3 )Ai(t), (5.9) 

where Mi, M2, M3 are constants depending on ||vo||, ||(/2o||j7 2 . Ml, Ms, A, X, e, but not on 

Proof. We note that in the following calculation only the lower-order uniform estimates (|3.1|) are 
used. The possible relaxation on the viscosity ^4 enable us to avoid using the L°°-norm of Vy, 
which was crucial in the proof of Lemma 13.61 In what follows, the generic constant C will only 
depend on ||v ||, H^oIIh 2 ) Mi, Ms, A, K, e. 

We revisit the terms on the right-hand side of (|3.18p . 



/ (v • V)v • Avdx < ||Av||||Vv|| L 3||v|| L 6 < C||Av||i||Vv||t 
< /z|ll Vv ll^ll Av H 2 + C M4^l|Vv|| 2 . 

h + ha + hb + hc<^- I (d i d i V J Ai) 2 ^ + C , ||Vv|| 2 3 ||Vd|| 2 6 ||d||i j0O , 

4 J Q 
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where 



l|Vv|i2 3 ||Vd||2 6 ||d||j> < C||Av|| 3 ||v|| 3 1| Vd|| 2 6 Hdlli,^ 

< C||Av||f ||v||3(||AV||5||A^||i + ||A^|| 2 )(||AV||||V^|| + ||V^|| 

< C||Av||i||v||3(||AQ||3 + l)(||e|| + l) 



< 
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+ Mj(i + lien s 



||Av|| 2 + C/V (||Q|| 2 + ||Vv|| 2 )||AQ|| 2 + C/i^UVvll 2 + C^\\Qf 



Next, 



hd < C||d||3 0O ||Ad||||Vv|| 2 4 <C'||V V 9||| 2 ||V V 9||| [1 ||VA^||||Vv|| 2 4 

< C(||VA^||I + l)||Av||l||Vv||5 < C(\\Q\\* + 1)(||AQ||3 + l)||Av||l||Vv|| 



< 



^ +M |(livv||i + ne|p- 

+C(^ 4 f +^ 4 3 )||Vv|| 2 . 



||Av|| 2 + C( M4 ^||Q||+^ 4 - 3 ||Vv|| 2 )||AQ|p 



h + h < ||Av||||vv|| L 3||d|| L oo||vd|| L6 < iiavI^iivv^iiv^ii^iiv^ii^iiv^iIh^ 

< C(||VAp||§ +l)||Av||l||Vv||i 
It is easy to see that 1% + I4 can be bounded just like l2d, because its order is lower. 



Av.(Qd)dx< ||Av||||Q||||V^|| LO c < ||Av||||Q||(||Q|| z + 1) 

T 3 

(g + ZxIIIQII 1 ) IIAvf + C^ + ^HQf. 



For the terms J±, J5 on the right-hand side of (|3,2ip . we have 

Ji < K||V^|| L oo||Av||||AQ|| +C||AQ|| || Vv|| L3 ||^|| W 2,e + C||AQ||||v|| L 6|M| w s,a 
< <7(||Q||* + l)||Av||||AQ|| < ^\\Avf + Oil Hi + IIQIDIIAQII 2 . 

J2 + J4 < C||VQ||(||V^|| 2 j o + l)(||Vv||||V V 9|| L - + ||v|| L 6||VV^|| L3 ) 

< q|VQ||||Vv||(||Q||i+l) 

< C(||AQ||3||Q||I + ||Q||)||Av||3||v||l||Q||! + C(||AQ||3||Q||I + ||Q||)||Vv|| 
XK 



< 



+ ||AQ|| 2 + / u 4 2 ||Q||||Av|| 2 + C(l + / u 4 



+ C||Vv|| 2 . 



XK, 



J 3 + J 5 < C||VQ|| 2 3(||V< 6 + 1) < G(||AQ||i||fi||5 + HQII 2 ) < — ||AQ|| 2 + C\\Q\\ 
Collecting all the estimates and using the Young inequality, we can obtain that 

4(l|Vv|| 2 + ||Q|| 2 ) + /ii f (didjVtDijfdx + Afis [ (d k ViD ki ) 2 dx 
at Jq Jo 



+ 



^-C/iKl + llVvf + IIQII 2 ) 



||Av| 



+ XK- C^(l + H 5 )(1 + ||Vv|| 2 + HQf) 
< C(l + /i 4 3 )(||Vv|| 2 + ||Q|| 2 ), 
which yields (|5.9p . 



l|AQ|p 



□ 
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Based on Lemma I5.2| one can prove the existence and uniqueness of global strong solutions 
(v, cp) to our system provided that the viscosity fj,^ is properly large. 

Theorem 5.3. Suppose n = 3. For any (vo,0o) € V x H 4 , if /Z4 > /x (vo,yo) ( c f- (|5.1ip ), 
problem (|1.5p ~ (|1.7p admits a unique global strong solution. 

Proof. The crucial step is to obtain a uniform bound of A\(t). Without loss of generality we 
assume that fj,4 > 1. Then we deduce from (|5,9p that 

d ~. 



dt 



Ai(t) + (y " Mi/xfAi^J ||Av|| 2 + (xK - 2M 2 /i 4 *Ai(t)j ||AQ|| 2 < 2M 3 A 1 (t). (5.10) 



ipr2| yields that Ai (r)dr < / + °° Ai (r)dr + 1 < max {2, £(0) + 1 =: M. If the viscosity 
/X4 satisfies the following relation 

ii 4 > /i 4 := max{l,K 2 }, with « := max ^2Mx, (Ai(0) + 2M 3 M + 2M). (5.11) 

then applying the classical method in [2D], we can argue as in |32| to obtain that 

y -Mi/j|Ai(t) > 0, Ai^-2M 2 ^4^Ai(t) > 0, V t > 0. 

The proof is complete. □ 

Finally, we study the long-time behavior of global solutions. 

Lemma 5.3. Let n = 3, the weak (or strong) solution (v,<p) to problem (|l,5p - (|1.7p has the 
following property: 

lim (||Vv(*)|| + ||Q(*)||)=0. (5.12) 

t— >+oo 

Proof. Since we are only concerning the behavior of (v, ip) for large time, due to the eventual 
regularity of weak solutions, we can reduce to the case of strong solutions by a finite shift of 
time. Then we can see that ||Vv(t)|| and ||Q(t)|| are uniformly bounded for t > 0. It follows 
from (|in$ that j|A(i) < C (similarly, from ([53p . we have ^Ai(i) < C). Recalling that 
A(t), Ai(t) € L 1 (0, +00) (cf. (|3.2p ). we arrive at the conclusion. □ 

Based on Lemma 15. 3| we are able to prove the convergence to equilibrium result in 3D. 
One can check the argument for 2D case in the previous section step by step. By applying 
corresponding Sobolev embedding theorems in 3-D, we can see that all calculations in Section 
4.2 are valid. Hence, the details are omitted here. 

Remark 5.2. Since the set of equilibria can form a continuum, the global solution obtained in 
Corollarv \5.1\ or in Theorem 1 5. 2\ will converge to an equilibrium (p^ which is not necessarily the 
original minimizer ip* . However, we can show that Efaoo) = E(<p*). To see this, we recall the 
definition ofw in the proof of Theorem \5. 6 A Actually we showed that the solution ip(t) will stay in 
the H 2 -neighborhood of ip* with radius less than (3, so does ip^. Then, we can apply Lemma\J^J 



with = tp* and ip = ipoo obtaining that \E((p 00 ')~ ^'( ( / 9 *)| 1 d — II — KA 2 ip 00 + V-f('S/ip 00 )\\A = 0. 
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